Kinetic accessibility of buried DNA sites in nucleosomes 
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Using a theoretical model for spontaneous partial DNA unwrapping from histones, we study the 
transient exposure of protein-binding DNA sites within nucleosomes. We focus on the functional 
dependence of the rates for site exposure and reburial on the site position, which is measurable 
experimentally and pertinent to gene regulation. We find the dependence to be roughly described 
by a random walker model. Close inspection reveals a surprising physical effect of flexibility-assisted 
barrier crossing, which we characterize within a toy model, the 'Semiflexible Brownian Rotor'. 



Although the DNA in eukaryotic cells is packaged into 
chromatin, its genetic information must be accessible to 
proteins for read out and processing The structural 
organization of chromatin is fairly well known: the funda- 
mental unit is a nucleosome core particle (NCP) consist- 
ing of about 150 base pairs (bp) of DNA wrapped in 1.7 
turns around a cylindrical histone octamer Q , and NCPs 
are regularly spaced along the DNA, which is further 
compactified into higher order structures. In contrast, 
the conformational dynamics of chromatin is poorly un- 
derstood. Recent experiments studied these dynamics on 
the level of individual NCPs using single-molecule force 
Q and fluorescence [1,0] techniques. The latter directly 
observed spontaneous conformational transitions where 
part of the DNA unwraps reversibly, allowing proteins to 
access DNA sites that are normally buried. This mode 
of access, driven by thermal fluctuations, is particularly 
important for passive DNA-binding proteins, e.g., tran- 
scription factors. Here, we study spontaneous DNA un- 
wrapping within a theoretical model, see Fig.[l](a). 

Consider a buried DNA site that is accessible only 
when a DNA segment of length L is unwrapped. How 
long is the typical dwell time Ta in the accessible state, 
i.e., the window of opportunity for protein binding? And 
what is the typical time ri for which it remains inac- 
cessible? Li et al. [3| measured Ta = 10 — 50 jms and 
Tiw250 ms for L^30 bp, while Tomschik et al. j5| found 
Ta = 100-200 ms and Ti = 2-5 s for L~60 bp. Taken to- 
gether, these results indicate a significant dependence on 
L in both time scales, which cannot be reconciled with 
an early theoretical study suggesting an all-or-none 
unwrapping mechanism where the nucleosome fluctuates 
between two conformations only. Instead, these results, 
as well as previous biochemical experiments f?], imply a 
multistep opening mechanism. 

In this Letter, we propose and characterize a theoreti- 
cal model for this multistep mechanism, similar in spirit 
to previous work on histone-DNA interactions which fo- 
cused mainly on static properties or the calculation of 
free energy barriers ^. Within our model, we clarify 
the physics that determines the L dependence of the time 
scales Ta and n. We find that the dependence of n can 
be interpreted with a simple random walker model, which 



may serve as a fitting model for future experiments that 
probe the time scales at different L values. In contrast, 
the L dependence of Ta reflects the intricate coupling be- 
tween the DNA polymer dynamics and the dynamics of 
breaking and reforming DNA-histone contacts. To ana- 
lyze the effect of this coupling, we introduce a toy model, 
the Semiflexible Brownian Rotor (SBR), see Fig. W[h). 
We identify a generic physical effect of flexibility-assisted 
barrier crossing, which may arise also in other contexts. 
It is marked by a characteristic plateau of the time scale 
at intermediate L. Biologically, the L dependence is rele- 
vant, because it creates a positioning effect for transcrip- 
tion factor binding sites relative to nucleosomes . We 
expect that the integration of single NCPs into nucleo- 
some arrays will alter the absolute time scales but not 
the basic physics of the DNA (un)wrapping process. 

Nucleosome model. — The NCP crystal structure Q 
shows that both the electrostatic and hydrogen bond in- 
teractions between the DNA and the histone complex 
are mainly localized to f4 contact points, about evenly 
spaced by 10.2 bp along a superhelical contour with ra- 
dius 4.2 nm and helical pitch 2.4 nm. Because we are 
interested only in the dynamics at a fixed (physiological) 
salt concentration, we combine the interactions at each 
of these points into a simple Morse potential [lH. The 
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FIG. 1: (a) Illustration of our nucleosome model. The DNA- 
histone interaction is localized at contact points attracting 
the red (dark) beads. The DNA is shown in the ground state 
as well as a conformation where the first contact is open, (b) 
Illustration of the Semifiexible Brownian Rotor (SBR) model. 
In this toy model, the tradeoff between bending energy and 
DNA-histone interaction in the nucleosome is mimicked by an 
angular potential V{ip), exerting a torque on the attachment 
angle of a semiflexible polymer at the origin. 
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DNA-histone interaction energy is then 

f/c = 7fcBT^(l-e-l'^'(")-'^"l/'')' 



(1) 



where c„ is the nth contact point on the superheh- 
cal contour, 7 is the depth, and p the width of the 
contact potential. A discrete bead-spring model with 
beads at positions models the DNA, and is 
the bead bound to contact n in the fully wrapped 
state. The beads are connected by a harmonic poten- 
tial Ub = eB^i{\Yi+i—Vi\ — af/2 with a typical bead 
separation a and a stiffness Eg set to SOO/csT/nm^. Be- 
low, we use 3 beads between contacts and at each end 
(about 2.5 bp/bead), unless stated otherwise. Increas- 
ing the discretization or Eg raises the computational ef- 
fort without affecting our results qualitatively. We ac- 
count for the bending rigidity of DNA by an energy 
Uh = e^Y^^ {I — cos 9i) with bending angle 9i at bead 
i and a bending stiffness eh adjusted such that the 
apparent persistence length matches the known £p « 
50 nm for DNA at physiological salt conditions. Fur- 
thermore, we incorporate the screened electrostatic self- 
repulsion of DNA through a Debye-Hiickel potential 



r,| with the 



Bjerrum length Is ~ 0.7 nm, a charge density t — 2 
charges/bp, and a screening length « 1 nm. We 
use a contact radius p = 0.5 nm in between the range of 
hydrogen bonds and electrostatic interactions and adjust 
the depth 7 of the Morse potential to match the bind- 
ing free energy [I4I of « 1.5A:RT_per contact estimated 
from biochemical experiments [7|,|9|. Taken together, the 
total energy is U — Us + Uh + Ubh + Uc- To study the 
dynamics of our model, we perform Brownian dynamics 
simulations with the overdamped Langevin Eqs. 



m = -MbVr.C/({r,})+r7,(i) 



(2) 



where pb is the bead mobility, and the absolute time 
scale is set by a? / pbhsT. The random forces r}^ satisfy 
{'n^{t) ■ r)j{t')) = 6pbkBT6i^j5it - t'). 

Unwrapping dynamics. — A suitable reaction coordi- 
nate for the opening of a single contact is the attach- 
ment angle see Fig. [Ija), which changes by Atp « 45° 
in this process. The equilibrium distribution p{'^) for 
the first contact is shown in Fig. [D^a) . Its bimodal form 
suggests to approximate a contact by a 2-state system, 
with rates kh, fcu for binding and unbinding, respectively. 
To test whether such a reduced description is sufficient, 
we initiate simulations in the fully wrapped state and 
determine the functionally relevant time scales, i.e., the 
average time Ti(n) until contact n opens to expose the 
nth DNA segment and the average time ra(n) until con- 
tact n recloses , • The results are shown in Fig. El^b) 
for n < 5 [3] ■ Within the reduced description of consec- 
utive 2-state contacts, ri(n) can be calculated as a mean 
first passage time jl6| for a ID biased random walker 
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FIG. 2: (a) Equilibrium distribution of the DNA angle de- 
fined in Fig. [TJa). The two peaks at = and tp ^ 45° 
correspond to the fully wrapped state and the state with con- 
tact 1 open, respectively, (b) Kinetics of DNA site exposure 
within our nucleosome model. The dwell time in the inac- 
cessible state (squares) increases roughly exponentially with 
the number of contacts that must open to render a DNA site 
accessible. The dashed line is a fit to Eq. (|3]). The circles 
show the average time the nth contact point remains open. 



with hopping rates fcu, fcb- The walker starts at site zero 
(reflecting boundary) and reaches site n after an average 
time 



Ti(n) 
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(3) 



Here, K — kh/k^ can be interpreted as the effective equi- 
librium binding constant per contact. The exponential 
increase of ri(n) is clear also from the equivalence of the 
biased random walk with a random walk against a free 
energy ramp. The excellent fit of ([3]) to the simulation 
data (dashed line) indicates that the reduced description 
is sufficient for the dwell times in the inaccessible state. 
In contrast, it proves insufficient for the dwell times in 
the accessible state, because ra(n) in Fig.^b) is clearly 
not constant as one would expect with a fixed binding 
rate kh- Thus, we find Ta(n) to be a more sensitive probe 
for the physics of spontaneous site exposure than ri(n). 

To probe the effect of the DNA length on the rewrap- 
ping kinetics, we vary the number of overhanging beads 
before contact 1 and plot Ta(l) as a function of the over- 
hang length L in Fig.[3l^a). Superimposed is the data of 
Fig. [2l^b) (bottom) with n converted to contour length. 
The good agreement of these dependencies indicates that 
Ta is determined by polymer dynamics. Indeed, we will 
now see that contact breaking and reformation of a ro- 
tating semiflexible polymer displays much richer physics 
than a simple ID barrier crossing process. 

Semiflexible Brownian Rotor. — The essential physics 
of contact formation in the nucleosome is captured by 
the toy model depicted in Fig.[Hb): A semiflexible poly- 
mer with contour length L and persistence length £p is 
attached to a point about which it can rotate in a plane. 
The attachment angle (p experiences a periodic poten- 
tial V{p) — Vb cos(27r iy9/A(/3), which creates preferred 
angles separated by potential barriers as in our nucle- 
osome model (there, the barrier for contact reformation 
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FIG. 3; (a) The dependence of the dwell time ra(n = 1) on 
the overhanging DNA length (diamonds) is compatible with 
Ta(n) when n is converted to contour length (gray circles). 
The dashed line indicates the diffusion limit (see main text 
for details), (b) The average barrier crossing time Tw (circles) 
for the SBR model of Fig. [ijb) . At small lengths, the barrier 
crossing time follows that of a stiff rod (indicated by the dot- 
ted line). Beyond a crossover length £c <^ £p, barrier crossing 
is much faster than for a stiff rod. For large lengths, Tw ap- 
proaches the diffusion limit, i.e., Tw of the free SBR (squares). 
With L < £p, free diffusion of the SBR (squares) is virtually 
indistinguishable from free diffusion of a rigid rod (dashed 
line). The crossover from the rodlike regime to the intermedi- 
ate regime is well described by the theoretical analysis (solid 
line), see main text. 



results from the DNA bending energy and the electro- 
static repulsion). The main difference is that the length 
of the rotating polymer is constant for this 'Semiflexible 
Brovifnian Rotor' (SBR), while it changes slightly when 
a contact breaks or reforms in the nucleosome. Also, we 
do not consider a directional bias in the SBR, because it 
is not essential for what follows. So far, barrier crossing 
of semiflexible polymers was studied only for situations 
where the entire polymer experiences an external poten- 
tial 17]. In the NCP, the potential acts only on the angle 
at the attachment point. 

To characterize the phenomenology of the SBR, we de- 
termine its barrier crossing rate 1 /r^ with Brownian dy- 
namics simulations of a discrete bead-spring model [iSj . 
The circles in Fig. ^h) show as a function of L/ip for 
Vq = bkBT. We observe that at very short lengths, Tw 
follows the stiff rod behavior [l^ indicated by 

the dotted line. However, above a certain length £c, there 
is a regime where is nearly insensitive to L, before it 
rises again. Hence, for lengths L > £c the semiflexible 
polymer crosses the barrier much faster than the stiff 
rod. What is the physical mechanism for this accelera- 
tion? One effect of a finite flexibility is a reduced mean 
end-to-end distance (due to the undulations in the con- 
tour), which in turn leads to a larger rotational mobility. 
However, with V{(p) — 0, the rotational diffusion time 
of a semiflexible polymer over an angle Atp (squares) is 
almost identical to that of a stiff rod (dashed line) when 



L < £p. Hence the acceleration is not a mobility effect. 
Note that the dashed line is also the diffusion limit for 
Tw, which induces a second crossover from a reaction to 
a diffusion controlled process. The equivalent diffusion 
limit is shown also in Fig.[3][a) (dashed line). It indicates 
that the Ta,{n) data for the nucleosome is indeed in the 
accelerated barrier crossing regime. 

Flexibility-assisted barrier crossing. — To understand 
the interplay between the polymer dynamics and the bar- 
rier crossing dynamics qualitatively, we recall the basic 
aspects of each: (i) A semiflexible polymer of length L 
relaxes its conformational degrees of freedom in a time 



Ly£p 



Conversely, within a given time r, a 



local bending deformation is "felt" only over a length 
£ ^ {£pTy/'^. (ii) The probability current over a bar- 
rier is proportional to the quasiequilibrium occupancy of 
the transition state and to the relaxation rate out 
of this state. Together, (i) and (ii) imply that £c is the 
length of the polymer segment that gets deformed during 
the relaxation process away from the potential peak. We 
estimate £c by noting that the attachment angle relaxes 
according to (p = —fi{£c)dV/d(p, where fi{£c) ^c'^ is 
the rotational mobility of the deformed segment. Hence, 
r - £3(^^/27r)Vt/o and with 4 ^ (^r)!/*, we find 



£r — C £rs 



"27 



(4) 



where C is a constant to be determined below. For 
lengths below ^c, the entire polymer is involved in the 
relaxation process, i.e., it behaves like a stiff rod. 

Quantitative theory for the crossover. — To render the 
above picture quantitative, we employ the Langer the- 
ory for multidimensional barrier crossing processes 21 1 . 
For the case at hand, one can show [2^ that the bar- 
rier crossing time simplifies to Tw = j!-e^^°/'^^^, where 
A_ is the eigenvalue associated with the unstable mode 
at the saddle point. We calculate A_ using the contin- 
uous wormlike chain model in the weakly bending ap- 
proximation [23| . At the transition state the chain is 
straight, e.g., along the x axis. We denote deviations 
from this configuration by y{x,t). The chain dynamics 
follows dty = — (fcsT£p/C) d^y with a friction coefficient 
C- With F = Vo{2tt / Aip)"^ denoting the curvature of the 
potential at the transition state, the torque on the at- 
tached polymer end is —T dxy\x=o- This torque must be 
balanced by a local bend resulting in the boundary con- 
dition kBT£pdxy\x=o — dxy\x=o- The other bound- 
ary conditions are y\x=o = dly\x=L = d^ulx^L = 0. 
We find a unique unstable mode with eigenvalue A_ = 
kBT£pa* /Ai^L* and a determined by 



Q;[sinh(Q!) — sin(a)] L 
cosh(Q;) + cos(a) +2 £c 



(5) 



where £c is as in ^ with C — s/T2. In the limit 
L £c, we find A_ = 3r/(L^ independent of the stiff- 
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FIG. 4: Dynamics at the barrier, (a) The unstable eigenmode 
for three different lengths. Polymers shorter than ic rotate 
without significant deformation, while long polymers form a 
bulge of size ~ Ic at the origin, (b) The prefactor of the 
Kramers time -?„ = 1/A_ as a function of the length. The 
prefactor increases as if L <g £c and is constant if L 3> 



ness, whereas in the opposite limit A_ = 3T/(£l indepen- 
dent of L. Fig. [4] shows (a) the unstable eigenmode for 
L/£c = {0.1,1,10} and (b) the crossover in the barrier 
crossing time. The eigenmode shape confirms our qual- 
itative picture: stiff and short polymers respond to the 
torque by rotating as a whole, whereas the torque shapes 
a bulge of size ~ £c in longer polymers. For a discrete 
polymer model, the same analysis can be performed, but 
the eigenvalue A_ must be computed numerically. The 
solid line in Fig. [3] shows the resulting barrier crossing 
time for the same discretization as used in the Brow- 
nian dynamics simulations of the SBR model. Indeed, 
the crossover from the rodlike to the flexibility-assisted 
barrier crossing is well described by this analysis. The 
deviations at larger L can be attributed to finite barrier 



corrections [24 1 



Discussion and outlook. — The experiments 0, IE] have 
shown that the functionally relevant time scales Ti and 
Ta depend on the position on the nucleosomal DNA. Our 
theoretical study suggests that these time scales addi- 
tionally depend on the total DNA length. The position 
dependence of ri should follow the random walker model 
(l3|), which is the minimal model for a gradual, multistep 
opening mechanism. However, we expect that the posi- 
tion dependence of Ta and the length dependence of both 
time scales will reflect the polymer dynamics of the DNA. 
Within our toy model, the Semiflexible Brownian Rotor, 
we find three physically distinct regimes for this length 
dependence, see Fig. [3ljb). The intermediate regime dis- 
plays a striking flexibility- assisted barrier crossing effect, 
the onset of which is marked by the new length scale 
£c of Eq. ([4]). It can be interpreted as the length over 
which the polymer contour is deformed as it passes over 
the potential barrier. Because £c is considerably smaller 
than the persistence length £p, we expect that the on- 
set of the intermediate regime will not be detectable 
in nucleosomes. However, nucleosomes should display 
the crossover from flexibility-assisted barrier crossing to 
diffusion- limited dynamics as shown in Fig. [3ja). All 
three regimes of Fig. could be probed in an experi- 
mental realization of the SBR model, e.g., with an actin 



filament as the rotating polymer. 
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